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Abstract 

This paper studies control problems of sampled data systems which are subject to random sample 
rate variations and delays. Due to the rapid growth of the use of computers more and more systems 
are controlled digitally. Complex systems such as space telerobotic systems require the integration 
of a number of sub-systems at different hierarchical levels. While many sub-systems may run on a 
single processor, some sub-systems require their own processor or processors. The sub-systems are 
integrated into functioning systems through communications. Comunications between processes sharing 
a single processor are also subject to random delays due to memory management and interrupt latency. 
Communications between processors involve random delays due to network access and to data collisions. 
Furthermore, all control processes involve delays due to causal factors in measuring devices and to signal 
processing. 

Traditionally, sampling rates are chosen to meet the worst case communication delay. Such a strategy 
is wasteful as the processors are then idle a great proportion of the time; sample rates are not as high 
as possible resulting in poor performance or in the over specification of control processors; there is the 
possibility of missing data no matter how low the sample rate is picked. 

Randomly sampled systems have been studied since later 1950*8, however, results on this subject are 
very limited and they are not applicable to practical systems. This paper studies asymptotical stability 
with probability one for randomly sampled multi-dimensional linear systems. A sufficient condition for 
the stability is obtained. This condition is so simple that it can be applied to practical systems. A design 
procedure is also shown. 


1 Introduction 

Many complex systems today involve the integration of a number of different subsystems at various hierar- 
chical levels. Examples of hierarchical subsystems are, for example, in the case of spacecraft: 

Level 1 - Assignment of systems to tasks; 

Level 2 - Assignment of subsystems to task systems, such as the shuttle manipulator, one of more cameras, 
an astronaut on EVA; 

Level 3 - Control of individual subsystems, cameras comprised of pan tilt, zoom, focus, feature tracking, 
exception warning; or control of machine tools comprised of spindle, table, tool changer, gauge; 

Level 4 - Control of elements, control of manipulator joints, end-effector force measurement, machine 
tool spindle drive, elevator motor drive, submarine plane control. 

These systems all comprise many components which may be ranked hierarchically. Many of the compo- 
nents are now computer controlled and are integrated by means of digital busses or networks. The integration 
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of these components into functioning feedback systems, camera - force-sensor - manipulator roll sensor, pitch 
sensor - main propulsion - planes, implies in addition to data communications, communication rates. 

In the case of communication networks, the data rates of point to point communication busses are well 
known. The rates at which a computer can respond to communication data interrupts requests add a variance 
to the data rates. In the case of shared networks, such as Ethernet, data collisions add considerable variance 
to the data rate frequently exceeding the data rate itself. However, these shared networks are very attractive 
from both the reliability and flexibility standpoints. 

Because of their flexibility in programming and speed in computing, digital computers are now regularly 
employed as integral components of dynamic feedback control systems. They are easily programmed to 
realize desired compensators. Due to the discrete nature of digital computers, variables in dynamic systems 
are sampled and quantized before sending to the computers. The well established discrete time system theory 
(e.g., [8]) provides methods to analyze the behavior of sampled data systems, based on the assumption that 
the sampling rates are fixed and the same, and the sampling operations on different channels of the systems 
are synchronized. If the sampling rates are fixed but different on different channels, known as multi-rate 
sampling, the system analyses are simple if the sampling rates have integral ratios [6, 10]. 

Due to random delays in measurement devices, signal processing, interrupt latency, priority scheduling, 
conditional branching, network communications, etc., sampling rates vary randomly in many systems, and 
the system performance could be expected to be improved if a theory supporting random sampling rates 
was used. Systems with random sampling processes are called randomly sampled systems. The behavior 
of a randomly sampled system is, presumably, related to the statistical properties of the random sampling 
processes as well as system parameters. Randomly sampled systems have been studied by Kalman [11], 
Leneman [16], Kushner and Tobias [15], Agniel and Jury [2], and others. One of the major motivations 
for studying randomly sampled systems in late 1950’s and early 1960’s was the introduction of digital 
computers in control systems. However as the speed of computers improved dramatically, time delays 
caused by computers became practically negligible in simple single processor controlled systems compared 
to other delays, and research on randomly sampled systems came to an end. Nowadays, development of 
computer controlled systems has reached beyond the stage of single processor control. Many subsystems are 
integrated into large systems. Furthermore, many complex dynamic systems impose demanding computation 
requirement. For example, computation time becomes a bottleneck in the implementation of dynamic control 
algorithms of multi-joint robot manipulators. Delay caused by computation and communication is no longer 
a negligible factor. 

Early researchers in the area of randomly sampled systems primarily considered stability conditions 
of the systems. Their work is briefly summarized below. Kalman carried out a comprehensive study of 
sampling systems [11]. He classified sampling into six categories: conventional sampling, nonsynchronous 
sampling, multiple-order sampling, multi-rate sampling, noninstantaneous sampling, and random sampling. 
For randomly sampled systems, Kalman showed that if the second moment of the output of an autonomous 
system is stable, the second moment of the output remains bounded when a bounded input is applied 
to the system. Based on his state space method [13], Kalman [12] also discussed the regulator problem 
and stability of a linear system described by independent random functions. This class of systems include 
randomly sampled systems. Thus the stability conditions obtained for this class of systems are applicable 
for randomly sampled systems. Kushner and Tobias [15] studied an autonomous linear system with linear 
and nonlinear feedback. Using a stochastic Lyapunov function, criteria for stability with probability one 
and s-th moment stability (s > 0) were obtained for scalar linear systems, and criteria for stability with 
probability one and second moment stability were obtained for multi-dimensional linear systems. Agniel 
and Jury [2] investigated asymptotic stability with probability one of a linear system with a saturating type 
nonlinear component. A computational procedure was provided to determine the largest stability sector of 
the nonlinearity for asymptotic stability with probability one. Using a stochastic Lyapunov function, Agniel 
and Jury in another paper [1] gave a condition for the asymptotic stability with probability one and the 
second moment asymptotic stability for single-input single-output multi-dimensional linear systems. They 
also showed that if an autonomous system exhibits asymptotic stability with probability one, the system is 
almost surely bounded input-bounded output. Leneman [16] studied a single-input single-output first order 
linear system with feedback. He derived the second moment of the output for the cases with and without 
input. The input is a stationary stochastic process independent of the sampling process. Consequently, a 
condition for the second moment stability was given. Assuming the independence of the sampling times and 
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the signals, Dannenberg and Melsa [7] took the expectation of a linear system equation, obtaining a system 
equation of expectation of the states and outputs. The first moment stability analysis is similar to that of 
deterministic sampled-data systems. An example of a spacecraft control problem was given, in which it is 
assumed that there is a probability of missing messages. The problem of random sampling of a random 
signal was studied by Bergen [4] and Leneman [17]. Their focus was on deriving expressions of the spectral 
density of a random signal after a random sampling. 

This paper studies the stability of randomly sampled systems in relation to the random sampling pro- 
cesses. Though Kalman [11] and Kushner [15] have obtained necessary and sufficient conditions for the 
stability in the second moment, it is not so easy to apply these conditions to practical systems. This paper 
studies asymptotic stability with probability one and gives a necessary and sufficient condition for one- 
dimensional systems and a sufficient condition for multi-dimensional systems. These conditions are easy to 
verify for given sampling distributions and are thus applicable to practical systems. 

In the next section, the asymptotical stability with probability one is defined. A sufficient condition 
is given for multi-dimensional linear time-invariant randomly sampled systems which is also necessary for 
one-dimensional systems. A design procedure todetermine feedback gains is obtained in Section 3. If we use 
a nonlinear compensator such as a computed torque controller for a robotic control system, then we would 
have a set of simple two-dimensional linear systems. In Section 4, the stability of such two-dimensional 
systems is considered and the design prcedure is shown for a Bernoulli distribution, a uniform distribution 
and a mixed uniform distribution. 


2 Stability 

Consider following linear time-invariant control system. 

x(t) = Ax(t) + Bu(t), (1) 

where x is an n-dimensional state vector, u an r-dimensional control vector, and A and B are n x n and 
n x r matrices, respectively. For this system, we apply a constant state feedback input 

u(t) - I<x(t k ), (2) 

from t = t k to t = <*+i( = t k + At), where I\ is an r x n matrix. Then a:(ijb+i) is given as follows. 

*(<*+i) = (*(A») + *(A k )I<)x(t k ), (3) 

where 

f Ak 

<I>(Ajfc) = exp(AAjfc), and ^(At) = / exp(Ar)drB. 

Jo 

Sampling interval A* is assumed to be subject to some probability distribution function F(A) or distribu- 
tion density function /(A) and A,- and Aj(i ^ j) are statistically independent of each other. For simplicity, 
we write Eq. ( 3) as follows 

** +1 = r(A*)**. (4) 

In this paper, we use the following matrix norm which is compatible with usual Euclid norm for vectors: 

l|r|| = Mr*r)} 1/2 , (5) 

where T* is the conjugate transformed matrix and fx(r) denotes the maximum eigenvalues of the matrix F. 
Note, however, that while the stability of the system (1) or (3) is invariant under a similarity transformation 
of the state variables, the matrix norm depends on the transformation, namely in general 

||r||#||r- 1 rr||. 

The stability of randomly sampled control system Eq. (1) is defined as follows. 


343 



Definition 1 (Stability) The randomly sampled control system Eq. (1) is asymptotically stable with 
probability one if 


Prob[ lim ||x*|| = 0] = 1 

it— MX) 

for any initial state Xq, where ||x|| is the Euclid norm of vector x. 

Now we define the following notation: 

E[u)] : Expectation of random variable w, 

V[w\ : Variance of random variable u> } 

and assume that 

£[{log(||r(A)||)} 2 j<oo. (6) 

Then a sufficient condition of the asymptotical stability is given in the next proposition. 

Proposition 1 (Sufficient Condition) Randomly sampled control system (1) is asymptotically stable with 
probability one if 

E = £[log(||T -1 r(A)T||)] < 0, (7) 


We also have 

Prob[\\T- l x k \\ < HT-^oll exp {k(E + e)}] > 1 - , (8) 

for any e > 0, where V = V [log(||T“ 1 T(A)T||)] - 

< proof > Assuming xq / 0 without loss of generality, from Eq. ( 4) we have 


k — i 

iog(||r- 1 **||/||r- 1 *o||) < j;iog(||r- 1 r(A,-)r||). 

i= 0 


Then the proposition is easily proved by the statistical independence of A*’s and Thebyshev’s inequality. 

< end of proof > 

We note that for one-dimensional systems the condition stated in the above proposition is necessary and 
sufficient for the aymsptotic stability with probability one [14]. If the sampling interval is constant, the 
condition in Prop. 1 is also necessary for the asymptotic stability of multi-dimensional systems. 

Now we define 


7 (A) = log(||T- 1 r(A)T||), and 


9(A) = 



(9) 


then we have the following proposition. 

Proposition 2 

i. If the sampling rate A is subject to a Bernoulli distribution where A = a with probability p and A = (3 
with probability q = 1 — p, then the system is asymptotically stable with probability one, if 


py(a) + qy(fi) < 0. 


ii. If the sampling rate A is subject to a uniform distribution U[a , 0\, then the system is asymptotically stable 

with probability one, if 

g(a) < g((3). 

iii. If the sampling rate A is subject to U[a, /?] with probability £ and to U[p,u] with probability 1 - e, then 

the system is asymptotically stable with probability one, if 


. 9(0) ~ 9(a) 
0- a 


+ ( 1-0 


g(Q - gQO 


< o. 


The proof is straightforward, so we omit it here. 
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3 Design Procedure 

Next we discuss a design procedure of a feedback gain K and a matrix T in the following. Now, assume that 
system 

x(t) = Ax(t) + Bu(t) (10) 

is controllable, then it is well known that the discretized system 

Zfc+i = $(Ajb)xjb + (11) 

is also controllable for almost all sampling interval A* [5]. Then we can assign poles {X iy i = 1, 2, . . . , n} to 
system (11) if poles {A*} are symmetric with respect to the real axis. Here, we apply Hikita s pole assignment 
algorithm[ 9 ] to the randomly sampled control systems. 

[Algorithm] 

step (i) For given {A,}, find r- dimensional vectors i = 1,2, . . . , n, which makes matrix 
T( A) = [v\ : • • ■ : v n ] non-singular. Vector v, } s are given as follows where <3> = $( A) and 
= *(A). 

• if A i is a real number, then 

v i = (*-\ i i n y 1 ni- ( 12 ) 

• if Aj and A*+i are conjugate complex numbers a, i jf3i, then 

Vi = Vuti — Vfeifi+ij and Vi+i = Vuti + V2*&+i, (13) 

where 

Vu = {(^-ai/n) 2 + /?, ? /n} _1 (^-«Jn)^, and V 2i = {($ - aj„) 2 + 0, I n ] ~ 1 /?»'*' • (14) 

step (ii) Feedback gain K is given as follows . 

I<( A)=-[6 :---:£n]T(A)- 1 . (15) 

step (iv) Check the stability using Proposition 1 or 2. If not stable , return step (i) and try 
another {Xi} and/or A. 

It is easy t 6 show that for this T(A) and A" (A), we have 

||r- 1 (A)r(A)T(A)|| = max{|A,|}. (16) 

Hence we can use matrices T( A) and /\ (A) to calculate 7 (A) and ^(A), In the next section, we use notations 
7 (A, A) and </(A, A) for 7 (A) and < 7 (A), respectively, to show the dependence of the functions on A clearly. 

4 Two Dimensional Systems 

In this section, we consider control of robot manipulators. We view a robot manipulator as a component 
of a large system, such as a space station. The robot controller communicates with the other components 
of the system to achieve cooperative actions. Communication between components is considered to have a 
longer delay than that within a component. We assume that robot controller has an inner feedback loop 
which compensates the nonlinearity of manipulator dynamics and operates independently of the other part 
of the system. The robot dynamic system together with the inner feedback loop becomes a linear system. It 
is feasible to treat the robot manipulator subsystem as a linear system when integrating and communicating 
with the other components. For example, if we use the nonlinear feedback controller developed in [3], we 
have n (=DOF of manipulator) decoupled two-dimensional linear systems 

*(!)=[ 0 {]*«)+[?]»«> <"> 
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Figure 1: function j(0 , 1) and g(0, 1) 



where x(f) = (e,(<),e,*(<)) is the error vector for the i-th component of outputs and u(f) is the corresponding 
input for this component of outputs. If the task is specified in joint space (the joint space control), the i-th 
component of output is simply the displacement of the i-th joint and the error vector is composed of the 
joint position error and joint velocity error. 

We now study the asymptotical stability of this system under the random sampling rate. The corre- 
sponding discrete time system is easily obtained for a sampling interval A as follows. 


zjt+i = 


1 A* 
0 1 


x k + 


Afc/2 1 


«*• 


(18) 


We apply the algorithm given above to this system directly. Then we have the following proposition. 


Proposition 3 (PD Controller) Assume that { A* ) = {Aj,A 2} where Ai ^ A2, then we have 
K( A) = ((Ai -f A2 — A 1 A2 — 1)/A 2 , (Ai + A2 + A1A2 — 2)/(2A)), 


and 

t(A, A) = 1), 

where 6 — A/A. 

The proof is obtained by direct calculation. This proposition implies that the function 7(A, A) is the 
same as the function 7 (6, 1) if we use K( A) = {k p /A? ) k v /A) instead of /f(l) = (k p ,k v ). Therefore we have 
<;(A,A) = A<7(#,1) for the same K(A). This fact is very useful to design the feedback gain. This will be 
shown by examples. 

Fig. 1 shows 7(0, 1) and g(0, 1) for = 0.4 and A 2 = 0.7, where we have 


K(l) = —(0.18, 0.81), and T( 1) = 


-0.759 -0.934 ' 

0.651 0.333 ’ 


and was used to make the norm of column vectors of T matrix be equal to one. 

Example 1 (Bernoulli Distribution) Let’s assume that the sampling interval is subject to Bernoulli dis- 
tribution, i.e. A = a with probability p and A = /? with probahhty q, where a < /), 0 < p < 1, and q = 1-p. 
The sufficient stability condition is given as follows . 

pj(a/A, 1) + tf7(/?/A, 1) < 0. (19) 

Note that if A > /?/ 1.96(= A*) then the system is asymptotically stable for any a because 7 (M) < 0 f or 
any 6 < 1.96. But we are generally interested in the smallest A because it gives us the fastest response. 
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Fig. 1 shows that the function y(6, 1) reaches the minimum value —0.417 at 8 — 1.35. Let 0* be the point 
which satisfies the following equation. 

7(0*, 1) = - x 0.417. 

9 

Then it is clear that A must be greater than A mtn (= 0/8*) for Eq. (19). 

A suitable value of A can be found from the range A min < A < A* by a trial- and- error method using 
Fig. 1 or Table 1 which gives pairs of {61,62} such that 7 ( 0 *, 1) = 7 ( 02 , 1)* 

(i) Calculate a — —(q/p) 7 (/?/A,l). 

(ii) Find {^ 1 , ^ 2 } such that 7 (^ 1 , 1) = 7 ( 02 , 1) < a using Fig. 1 or Table 1. 

(iii) Check 6\ < a/ A < 02- U so, calculate K{A). If not so , go back to step (i) with another A. 

For example, if a = 10 msec , 0 — 30 msec, and p — 0.75, then 0* is about 3.64 and A min = 8.24 msec , 
while A* ~ 15.3 msec. If we select A = 11 msec then ^(0/ A,l) = —0.278 and a/A = 0.91. Therefore 
we can try the 6-th row of Table. 1 , and we have 0 1 = 0.84 < 0.91 < 02 = 1.68. Hence the system ts 
asymptotically stable for K = —(1488,73.64). 

Example 2 (Uniform Distribution) Now assume that A is subject to a uniform distribution U[q,/ 3]. 
The sufficient condition of the asymptotical stability with probability one is given as follows: 

g(a/A, 1) > g{0/ A, 1). 

The function g(8 , 1) has its minimum value at 0 = 1.96. Now we define A* = /?/ 19.6 and A min = 0/2.89. 
If A > A*, then the above sufficient condition is satisfied for any a. Therefore the system ts asymptotically 
stable if A > A*. On the other hand , if A < A m i„, then the above condition is not satisfied for any ot. 

Table 1 also gives pairs of {63,64} and the ratio 63/84 such that </( 03 , 1) = #( 04 , 1 ). If there is a pair 
{ 03 , 04 } such that a/0 > 83/84, then the system is asymptotically stable for the K(A) where A — a/ 03 . 
Therefor we can determine A easily using this table as follows: 

(i) Calculate a = a/0. 

(ii) Find a pair { 03 , 04 } in the Table 1 such that a > 63/84. 

(ii) Calculate A = a/83 an d ( A). 

Now assume that a = 10 msec and 0 = 30 msec, then we have A* = 15.3 msec, A mm = 10.38 msec, and 
a/0 = 1/3 > 0.273 in the Table 1. Therefore we can use a/A = 0.75 and A = 13.33msec. Hence the system 
is asymptotically stable with K = —(1065,62.31) if 0 < 36.7 msec. Table 2 shows the IAE (Integration of 
Absolute value of the Error) for fifty random streams with the initial condition z(0) = (1.0, 0) T . The table 
shows that when > 40 msec, the STD (STanderd Deviation) and the maximum values of IAE for the 
velocity error e.i(t) become very large compared to the cases where 0 < 35 msec. This means that the system 
is still stable but there is a large vibration in the response for A > 40 msec. It is interesting since A selected 
above assures the asymptotically stability for 0 < 36.7 msec. 

Example 3 (Mixed Uniform Distribution) Next we assume that A is subject to a uniform distribution 
U[a,0] with probability e and toU\}i, v\ with probability 1 - e. The sufficient condition is given as follows: 

E _ g(v/A,l) - gpi/A, 1 ) < Q 

0/A-a/A v/A-p/A 

Though the selection of A becomes a little difficult , we can use the following procedure to estimate an 
appropriate A: 

(i) Define a = (a + /?)/2.0, 0 = {p + ^)/2.0, p = e, and q - 1 - p. 

(ii) Determine A using the procedure in Exam. 1 for a = d and 0 — 0. 

(iii) Check the condition. If satisfied, calculate K( A). If not, try another value for A. 
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Figure 2: Simulations for Bernoulli Distribution, Uniform Distribution and Mixed Uniform Distribution 



time (in seconds) 


Table 1: 6 \, 9 2 , 63 , and #4 


7(*i,l) = T(*a,l) 

9(6 3,1) = 

--9(64,1) 

2l*i 11 

9 1 

62 

9(6 , 1) 

#3 

O 4 

03/04 

0.00 

0.00 

1.96 

0.000 

0.00 

2.88 

0.000 

-0.05 

0.18 

1.92 

-0.009 

0.25 

2.87 

0.087 

-0.10 

0.33 

1.89 

-0.039 

0.50 

2.84 

0.176 

-0.15 

0.46 

1.83 

-0.094 

0.75 

2.78 

0.270 

-0.20 

0.58 

1.79 

-0.173 

1.00 

2.69 

0.372 

-0.25 

0.71 

1.73 

-0.270 

1.25 

2.56 

0.488 

-0.30 

0.84 

1.68 

-0.373 

1.50 

2.39 

0.628 

-0.35 

0.98 

1.60 

-0.456 

1.75 

2.17 

0.806 

-0.40 

1.18 

1.48 

-0.467 

1.80 

2.12 

0.849 


Now assume that A is subject io ZV[5 msec, 15 msec] with probability e — 0.75 and to U[20 msec, 40 msec] 
with probability 0.25. Then we have a = 10 msec, /? = 30msec, p = 0.75, and q = 0.25. 7/we use A = 11msec 
from the result of Exam. 1, then we have E = —0.04 < 0. Therefore the system is asymptotically stable for 
the same K = -(1488,73.64). 

Fig. 2 shows the simulations of x(t) for three cases discussed above where x(0) = (1.0, 0) T . 

It is easily shown that even if we use a PID controller 

z k+ 1 = zjt + [ 1 : 0]xjt, and u k = + /f 2 xt, (20) 

or a PD controller with one step delay 

«* = /'f($(A)x*_i + ¥(A)u t _, ), (21) 

instead of the PD controller given in Prop. 3, we have the similar proposition . Therefore we can determine 
A easily. 
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5 Conclusions 


In this paper, the stability of randomly sampled linear control systems was discussed and the following results 
were obtained. 

1. A sufficient condition for the asymptotical stability in a norm with probability one was obtained for 
multi-dimensional systems. 

2. For a simple two-dimensional system with PD controllers, a design procedure was shown which was 
easily applicable to systems with PID controllers or PD controllers with one step delay. 

The results given in this paper are also easily applicable to the robotic control systems where computed 
torque controllers or PD controllers with a feedforward term are used at the random sampling rate. The 
results will be shown in the near future [14]. 


Table 2: IAE for U[10 msec,/? msec] 



MEAN 

STD 

MAX 

p 

mm 

BHH 

iK&a 

KSfiBl 

e,(<) 

EB2H 

25 

0.0531 

0.9999 

0.0022 

0.0011 

0.0572 

1.0020 

113811 

0.0511 

1.0039 

0.0043 

0.0200 

0.0561 

1.1310 

KH9 

0.0498 

1.0198 

0.0060 

0.0515 

0.0589 

1.2370 

ran 

0.0509 

1.2418 

0.0077 

0.4305 


mm 

lai 

0.0709 

2.6492 

0.0638 

4.4592 

■flEKfjfl 

wmm 
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